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Abstract 



We consider a model of three electrons and one hole confined in a 

two-dimensional (2D) plane, interacting with one another through Coulomb 

forces. Using a Ritz variational method we find an upper bound of ~ 

— 0.0112me 4 /87r 2 e 2 7l 2 for the ground-state energy of such a system when the 

particles are near one another. The possible connections of such a complex 

to other fields of physics are discussed. 
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I. INTRODUCTION 



The research on trions in 2D systems has shown many interesting results in recent 
yearsEHi. In a broader sense, the trion problem, not limited to 2D, is linked to the existence 
and the stability of negative hydrogen ion H~ which has not only important implication for 
atomic physics and astrophysics!, but also vital application to other fields such as physical 
chemistryil, accelerator engineeringB, and controlled fusion!. 

In a trion, or charged exciton, two electrons are bound to one hole despite the strong 
Coulomb repulsion between electrons. (Likewise two holes can be bound to one electron 
as well but the discussion that follows is the same if "hole" is switched with "electron.") 
It is quite a marvel that the two electrons would not fall apart from the strong Coulomb 
repulsion between them. More surprising is that the ground-state energy of such a complex 
is lower than that of an exciton, at least in the limit of heavy hole (i.e., the mass of the hole 
is much larger than that of electrons, which effectively results into the H~ problem). 

First let us try to understand in terms of classical mechanics why trion forms despite the 
strong Coulomb repulsion. 
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FIG. 1. Two electrons and a hole arrange themselves to form the bound state of trion in a 
classical picture of electrodynamics. Arrows indicate forces. Electrons el and e2 are shown as the 
shaded circles. 

Suppose that two electrons are on two sides of one hole as shown in Fig. 1. The 
electrostatic forces felt by electron 1 (el) are from the hole and electron 2 (e2). Because 
the forces are governed by the inverse square Coulomb's law, the repulsion exerted on el by 
e2 is only 1/4 that of the attraction given by the hole. Similar argument applies to e2 and 
therefore these two electrons tend to come together to form a bound state in spite of the 
strong repulsion among like charges. In order for electrons not to fall into the hole we can 
assume in a classical fashion that the electrons rotate around the hole in such a way that 
the inward electrostatic force is balanced by the centrifugal force. In quantum mechanics 
such tendency of electron's coming together will be stopped by the penalty in kinetic energy 
which increases as a -2 , where a is roughly the average distance between the hole and one of 
the electrons. Note that the gain in Coulomb energy varies at a slower pace approximately 
as a -1 so that it would be hard for the electrons to come near the hole without limit, if the 
possible electron-hole recombination is not taken into account. (The trion will not be stable 
if the recombination is spontaneous. We will discuss this problem further in Sec. [VI] .) 

Therefore the bound state of trion is not quite a surprise after all. It basically results from 
the inverse square Coulomb's law and a favorable geometric configuration of two electrons 
and one hole. 
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II. THREE ELECTRONS AND ONE HOLE 



Now we want to take one step further to see if a similar physical picture generalizes to 
the case of three electrons and one hole. In the heavy hole limit, this problem is in effect that 
for the doubly charged negative hydrogen ion H . There is a much celebrated theorem!!^ 
by E. H. Lieb saying that H is not stable against dissociation into H~ and e~ in vacuum. 
Here vacuum means that there is not any other charged particle nearby at finite distance 
except for the particles forming the H ion. We need a few more words on this theorem. In 
Lieb's original papers, a bound system is defined as an eigenstate of the Hamiltonian. This 
is a strong condition, which leaves room for metastable states for H . The existence of 
possible metastable states for H is currently a highly debated issu Lieb's theorem 

either does not rule out any possible H state whose total energy is negative. Here the 
energy zero is taken when particles are infinitely separated and thus with arbitrarily small 
kinetic and potential energy. 

In this paper we study the bound on the ground-state energy of three electrons and one 
hole at close distance in 2D. Instead of having a heavy positive charge, we assume that the 
masses of the four are all the same. Our result shows strictly that the ground-state energy 
is negative when the four particles are relatively localized one to another. 
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FIG. 2. Three electrons and a hole arrange themselves to form a "bound" state in a classical 
picture. Arrows indicate forces. The dotted equilateral triangle is a guide to the eye, showing three 
electrons at equal distance. The dash-dotted circle traces the orbit of electrons. 

Before starting quantum mechanical calculation, we want to show in a classical picture 
why such a state of negative energy might be possible. In Fig. 2, three electrons are at 
the vertices of an equilateral triangle, with one hole in the center. According to the simple 
knowledge of electrostatics, the force on the central hole is while those on the electrons are 
4v%~~^ inwards to the center, where a is the distance between the hole and one electron. 
Not to make themselves fall into the hole the electrons may rotate around the hole in such a 
way that the inward electrostatic force is balanced by the centrifugal force, just as in a trion. 
Intuitively a four-body system with a net charge of — 2e may not be stable because of the 
strong Coulomb repulsion. However under the special geometry as shown, the peculiarity of 
Coulomb's law (an inverse square law) actually favors a bound state of 3 electrons and 1 hole 
(3E+1H), at least within the framework of classical electrodynamics. Thus these electrons 
tend to move to the hole, giving rise to an effective binding among them. 

Note that the total charge of such a complex is — 2e which just reminds us of the Cooper 
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pair that brings in superconductivity, because the complex of four fermions (three electrons 
and one hole) is now effectively a charged boson. 

The stability of such complex against electron-hole recombination or annihilation ac- 
tually depends on the energy band structure of the host conductors. While electrons and 
holes recombine spontaneously in semiconductors, they coexist peacefully inside semimetals 
without destroying one another. The stability is governed by semi-metallic energy band 



structures. We will deliberate on this later in Sec. [VI. 

For the purpose of convenience, let us introduce a new term, tetron, to refer to the 3E+1H 
complex, in a sense similar to the term trion. We will use tetron and 3E+1H interchangeably 
thereafter. 

Our strategy for this paper is as follows. We will first mention relevant 2D quantum 



mechanical results for the Coulomb potential in Sec. |TTT| . Then as a warm-up exercise we 
try to find out the upper bound on the ground-state energy of a 2D trion in the heavy hole 
approximation, which is given in Sec. [TV]. Based on techniques developed for trions, we try 
a few fully antisymmetrized wavefunctions in Sec. [V] to show that states of negative energy 
indeed happen to tetron. 

III. QUANTUM MECHANICS OF AN ELECTRON IN 2D COULOMB 

POTENTIAL 

This problem has been worked out numerous times but we just repeat some relevant 
results here for quick reference. 

The Hamiltonian for an electron in a 2D Coulomb potential reads 

H = -—V 2 --^—, (1) 

2m Aner 

where e is the permittivity, e the elementary charge and m the mass of the electron. 

We try to solve the problem in a polar coordinate, where the gradient operator is given 

by 
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^ * dr ^ r 89' ^ ^ 



with the corresponding Laplacian 



V 2 -i — (r— ) + — — (3) 

r dr dr r 2 dO 2 

The Schrddinger equation reads 

h 2 e 2 

{-—V 2 - — = m. (4) 

2m Aner 

We want to simplify Eq. (ffl) by making it dimensionless. Note that h 2 /2m has a dimen- 
sion of (energy) (length) 2 while e 2 /47re is in (energy) (length). We select the units (or the 
energy and length scales) so that both h 2 /2m and e 2 /4ire are 1. Under such a convention 
length will be measured in the unit 2^eh 2 /me 2 while energy is measured in me 4 /&ir 2 e 2 h 2 . 
This unit system is summarized in Table |. 
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TABLES 



TABLE I. The units chosen to make the Schrodinger equation dimensionless 



Dimension 


(energy) (length) 2 


(energy) (length) 


(length) (energy) 


(wavenumber) 


Ordinary unit 


h 2 

2 m 


47T<E 


2iteh 2 me 4 
me 2 8ir 2 € 2 h 2 


me 2 
2neh 2 


Dimensionless unit 


1 


1 


1 1 


1 



Now the Schrodinger equation takes a simple form, 

{-V 2 - = E^. (5) 
Eq. (^j) has the following solutions as well known, 

*„ m (r, 9) = H^EPiggg) , m = o, 1, 2, 3, ■ ■ ■ , (6) 

V Z7T 



where w nm (r) is given by 

(2Kr) m+1 / 2 exp(-Kr) K(n + m- 1/2)! 1/2 
Wnm(r) ~ (2^)1 [ n(n-m-l/2)! J 

X!Fi(-n + m+ ^;2m + l;2Kr). (7) 

Here xF\ is the confluent hypergeometric function, and 

113 1 

K — — , n m — 0, 1, n . (8) 

2n 2' 2' ' ' ' ' 2 v ; 

The energy E nm of eigenstate \l/ nm is given by 

For the purpose of convenience, we list the first four eigenstates as follows, 

*i,o( r » fl ) = \/-exp(-r), 
2 V vr 



/ 4 r 

*§,±iM) = \/^T^ rexp (~3) exp ( ± ^)- 



(10) 



s 



IV. UPPER BOUND ON GROUND-STATE ENERGY OF A 2D TRION 



We take this problem as a warm-up exercise. This is a very illustrative example showing 
how to balance between the Coulomb and kinetic energy. It also gives us some intuition in 
choosing good trial wavefunctions and on how to deal with two-body Coulomb integrals. The 
basic line of reasoning follows that of H.A. Bethel, E.A. HylleraasB and others0. Although 
there have been recent researches^ in the same direction, our methods in dealing with the 
integrals are somehow different and hopefully easier to understand. 

For the purpose of simplicity we take the mass of the hole as infinity. The Schrodinger 
for this problem is now clearly 

{-Vl-V 2 2 - --- + — }^ = E^. (11) 
ri r 2 ri 2 

where and V 2 are respectively short-hands for 

1 d , d . 1 d 2 

ri or i or i r{ o6( 

and 

Id d Id 2 

— ^— (r 2 ^— ) + ^^2- (13) 

r 2 or 2 or 2 r 2 cw 2 

Here subscripts refer to electrons such that 7*1, 9\ and r 2 , 82 are the coordinates of el and 
e2 respectively. 

We choose the following trial wavefunction 

* = exp[-A(r a + r 2 )][l + fc(r a - r 2 ) 2 ], (14) 

where A and k are parameters to be determined later by variational principle. The first term 
in this wavefunction suggests that el and e2 are both in Is states. Here Is state are not 
exactly the given in Eqs. ( |T0D but rather that with a screening constant A. The physical 
origin of A is clear: the Coulomb field of the hole felt by one electron is reduced (screened) 
by the presence of the other. For this reason we also require that A be positive and smaller 
than 1. The existence of the second term is to reduce the Coulomb interaction between el 
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and e2. Note that this term is small when r\ and r 2 are close which is another way of saying 
that electrons 1 and 2 tend to avoid each other. So why not just keep the second term only? 
The reason is that while this term is good for reducing the Coulomb energy it causes some 
penalty in kinetic energy when compared with the first one. Thus the viable solution is to 
combine the merits of the two terms, which results in the trial wavefunction flT3|). 

To facilitate our calculation we want better management on the trial wavefunction. Note 
that functions like exp[— X(r\ + r 2 )] are not normalized and would be a headache for the 
algebra and calculus to come. The solution is to generalize the eigenstates of Eqs. fllP)) to 
accommodate the screening effect. The first four solutions now read 

| Is) = W^Aexp(-Ar), (15) 



i 2s > = fe |A -i A2r|exp( -T>' < 16 > 



IV) = y^AVexpt-^expfiiS). (17) 

where we have used notations Is, 2s, and 2p ± in place of the cumbersome The meaning of 
these symbols should be self-evident given a background knowledge of hydrogen eigenstates 
in three-dimension (3D). It is not hard to verify that these wavefunctions are now normalized. 
The trial wavefunction is then instead written as 

# = [l + k(r 1 -r 2 )*]\ls 1 ls 2 ), (18) 

where |lsils 2 ) = |lsi)|ls 2 ) stands for 



, ^Aexp(-Ar 1 )y^Aexp(-Ar 2 ), (19) 
and |lsi) is the short-hand for 

J^Xexpi-Xn). (20) 

Similarly we define |ls 2 ). 

This trial wave function is symmetric in the spatial coordinates and in order to make 
the total wavefunction antisymmetric we only need to multiply ^ by a spin singlet state 
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«^i2 = ^75 [| T1I2) — I I1T2)] where the subscripts are electron indices with the superscript a 
indicating antisymmetry. Note that the spin wavefunction in this case will not affect the 
calculation of either (^/S^ 2 \^S^ 2 ) = or (^S^ 2 \H\^Sf 2 ) = (^\H\^) thus we can drop it 



for now. 


We also tabulate some 


often used integrals in Table [TT] to facilitate our 


calculation. 




TABLE II. 


Some useful integrals involving Is) state 




Integral 


<ls|r _:L |ls) 


(Is r Is) (Is r 2 Is) (ls|r 3 Is) 


(ls|r 4 |ls) 


Value 


2A 


A- 1 |A- 2 3A- 3 


15 \-4 
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The Ritz variational principle states 

EW = i (W = min -' (21) 

where 

= (*| - V? - V* 1 1*)- (22) 

r 1 r 2 r 12 

First we need to calculate (^l^), 

(tf = (lsils 2 |{l + 2k{r\ + r 2 2 - 2r x r 2 ) 

+A; 2 [r^ + r 2 + 6r 2 r 2 - 4rfr 2 - 4rirf]}|lsils 2 ) 
= 1 + 2A;A~ 2 + 4.5A; 2 A- 4 . (23) 

The first term on the right-hand-side (RHS) of Eq. ( p2;) is 

- V 2 |^) = (Vi^|Vi*), (24) 

where |Vi^) is to be understood as Vi|\l/), while (Vi^|, as a notation, is the Hermitian 
conjugate of |Vi\I/). This identity can be proved in rectangular Cartesian coordinate using 
integration by parts when \I/ is zero at infinity which our trial wavefunction satisfies. 
Because 

V 1 |*) = {f 1 — + 0!- — }|1 S1 1S 2 > 

or i r\ oV\ 
= {-A[l + k{r\ + r 2 - 2r x r 2 )] + 2k(n - r 2 )}f 1 |ls 1 ls 2 ) 

= -A*fi + 2k(n -r 2 )fi|lsils 2 ), (25) 

we have 

(Vitf |Vi# ) = A 2 (tf |tf ) + 4£; 2 (3A~ 2 - 2A~ 2 ) 

= A 2 + 2k + 8.5k 2 \- 2 . (26) 

Since \I/ is symmetric in r± and r 2 it should not be a surprise that the second term on 
RHS of Eq. H, (tf | - VH*> is equal to (tf | - V 2 |^). 
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The third term in {V\H\V) [i.e., Eq. fl22|)] is 

- — |tf> = -(lsils 2 |— {1 + 2k(r 2 i + A~ 2r x r 2 ) 

+k 2 [r{ + r A 2 + Qr\r 2 2 - Ar\r 2 - Ar 1 rl]}\ls l ls 2 ) 
= -(2X + Ak\- 1 + 9k 2 \- 3 ) (27) 

to which the fourth term — is equal because of the symmetry in n and r 2 . 

The last term of (^\H\^) is a two-body Coulomb integral involving elliptic integrals, 
and 

(tf I — = (lsils 2 |— {1 + 2k{r\ + r 2 2 - 2r x r 2 ) 
n 2 r 12 

+k 2 [r{ + r 2 + <or\r 2 2 - Ar\r 2 - 4rir|]}|lsils 2 ) 
= (lsils 2 |— {1 + 2k{2r\ - 2nr 2 ) 

+k 2 [2r\ + 6rfr 2 2 - 8rfr 2 } } \ l Sl ls 2 ) , (28) 

where we have used identities 

(l Sl l S2 |l!|l Sl l S2 ) = (l Sl l S2 |lL|i Sl i S2 ) 
ri2 ru 

(l Sl l S2 |Il|l Sl l S2 ) = (l ai ls 2 |lL|i Sl i a2 ) 
ri2 r 12 

3 3 
(l Sl l S2 |^|lSlls 2 ) = (lSil5 2 |^|l Sl ls 2 >, 

ri2 

(29) 

which are easy to understand because of the symmetry of these integrals in T\ and r 2 . 

As an example, we explicitly show in the following how to work out the first term on 
RHS of Eq. d|). 

(lsils 2 | — |lsils 2 ) 

9 roc roo r2ir r2ir i 



•/ poo poo rzn rzn 

= (-) 2 A 4 / ndn r 2 dr 2 ex P [-2X(r 1 + r 2 )] / / rffl 2 . , (30) 

7T JO JO JO JO J r 2 + r 2 _ 2rir2 cos Q u 

where Q\ 2 = 6i — Q 2 . 
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We want to change the integration variables in the integral 

2vr r2n 1 

de x / dd 2 

>] J° - 'r\ + r\ — 2rir 2 cos #12 



(31) 



from 9\ and 62 to 612 and 62. To understand how this is done, we refer to Fig. 3. Shown in 
Fig. 3 is a plane formed by 6\ and 82- The integration limits of 6\ and 62 in Eq. [31] or the 



D (0, 2%) 



to the square ABCD. 






E (27C.47C) 




C (271,271) 






e 2 

1 ► 
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B (27t,0) 



FIG. 3. Change of variable from Q\ to #12 in integral (31). 
Because the integrand in integral (|3l|) 

1 (32) 
lr\ + r\ — 2r\t2 cos Q\2 

is periodic in 9\ with period 27r, the integration over triangle ABC is the same as that over 
triangle DCE. Note that DCE is the translate of ABC by 27r in the Q\ direction. Thus 
we can now change the domain of integration from the square ABCD to the parallelogram 
ACED. Since the segments AC and DE respectively satisfy the relation Q\ — 62 = 612 = 
and 61 — 62 = 6\2 = 27r, if we transform integral variables to 612 and 62 along with the change 
in integration domain, we will have 
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2tt i-2-k 

dO-L / d0 2 
Jo 







1 



r\ + r| — 2rxr 2 cos # 
1 



12 



0(01, 2 



2tt /-27r 

n ./() wrf + r| — 2rir 2 cos 



where I «/a 1 f 2 1 1 is the absolute of Jacobian determinant 

1 9(02,012) 1 



d(6i 


62) 




001 

802 


80! 
8012 




1 


1 


d(e 2 , 


012) 




802 
802 


802 
8012 




1 






which is identically 1 in its absolute value. 
Thus integral (|31|) now satisfies 

r2ir r 2ir 



rZIT rZIT 

I dd 2 / d9 12 
Jo 



1 



r\ + r| — 2rir2 cos #i 2 
1 



2- / d9 12 

/(l - 'rf + r\ — 2rir 2 cos 6 
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Inserting Eq. ( |35"1) into Eq. (|30l) , we get 



(lsils 2 | — |lsils 2 ) 
ri 2 



7T 



/■oo roc r-lix 

A 4 / drx / dr 2 exp[2X(r 1 + r 2 )]r 1 r 2 / d6 12 
Jo Jo Jo 



Vf + r 2 — 2rir2 cos # 
Eq. ( j36| ) can be further simplified by the following change of variables, 
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7™i = r, 
r 2 = ur, 



with corresponding Jacobian 

9(ri,r 2 ; 



<9(r, w) 



*9ri 
dr 


dri 
du 




1 





8r2 
dr 


8r 2 
du 




U 


r 



This transformation changes RHS of Eq. (^) to 

1 



—A 4 / MG^i / dQ\2 , 
7T JO JO V 1 + M — 2uCOS^i2 



/ <irr 2 exp[2A(l + u)r] 
Jo 
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which can be first integrated with respect to r to give 

8 ,„ . r°° r 2 „ , ,„ i 



X roc '/ fZTT 

-A 4 2! / [— -fudu / d9 12 - (40) 

Ti Jo [ A(l + u) J Jo VI + « 2 -2ucos^i 2 V ; 



where 



2- u 4:K( ~ Au 



d9 12 , = " " vl - 2 " + " 2; (41) 

o ^l + u 2 -2ucos6 12 VI - 2m + u 2 V ; 

after integration over #i 2 , with K the complete elliptic integral of the first kind. Integral 
(pfl]|) is therefore reduced to 



-Au 



—A / ] 3 ; Kl - 2u+u2 Ldu (42) 

vr 7o L l + u J Vl - 2m + it 



which can be integrated numerically to give the final answer 



(lsils 2 | — |lsils 2 ) = 1.1780972A. (43) 

r 12 



Similarly we work out other terms in Eq. (^ 



2?" 7" 

2fc(lsils 2 |— — |lsils 2 ) = l.OSOSSSlfcA -1 , (44) 



and 



ri2 



2 2r| + 6r 2 r| - 8r 3 r 2 h 
fc (lsils 2 lsils 2 ) 

= 1.5186410A; 2 A" 3 . (45) 



Substituting Eqs. (^3|) ,(^i|) and (|45|) into (p8|) and combining terms from (|26|) , (|27|) , we 



find the expectation value of the Hamiltonian 
($\H\V) 

= 2A 2 + 4k + 17k 2 \- 2 

-2.8219028A - 6.9691649/cA" 1 - 16.4813590A; 2 A" 3 , (46) 

where the first three terms are from kinetic energy and the last three are Coulomb ones. 
We want to minimize the energy 
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E 



(V\H\V) 



2A 2 + 4k + 17fc 2 A" 2 - 2.8219028A - 6.9691649A;A- 1 - 16.4813590£; 2 A- 3 

l + 2A;A- 2 + 4.5PA~ 4 



(47) 



with respect to A and k, which has numerical results as follows, 



E 



-1.08147; k = 0.111270; A = 0.693001. 



(48) 



This value, as an upper bound of the ground-state energy for our specific trion, is lower 
than that of an exciton with heavy hole, which is —1. The trion is therefore more stable 
than exciton and it forms spontaneously given one exciton and an extra electron. 

Our calculation shows that the trion binding energy in two-dimension > 0.08147 2D 
Rydberg, larger than that of 3D one which is ~ 0.055 3D Rydberg. This is not a surprise 
because it is well known that lower dimensionality favors the formation of bound states. As 
a result, pure excitons in 2D have a binding energy 4 times as big as that in 3D. 

We are so far confined to heavy hole limit. It is interesting to see quite some reportMBil 
showing trion states more stable than those of exciton regardless of the electron-hole mass 
ratio although the author of this paper have not obtained such results using a slightly 
different but more rigorous Hamiltonian. Not reaching a low energy after a few trials in 
variational methods does not rule out its actual existence and this problem will not be 
further probed here because it is off our main topic for this paper. Such issue aside, we have 
been familiar with the techniques needed for further development and let us now investigate 
the possible negative-energy state of tetron or 3E+1H. 



V. THE VARIATIONAL TREATMENT FOR TETRON 



A. The Hamiltonian 



We assume that all particles have the same mass m. The full Hamiltonian is 



H 



2m 



(V^ + V 2 + V 2 + V^) 
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where subscripts 1, 2, and 3 are indices to the three electrons while designates the hole. 
r 12 is the distance between particles 1 and 2. The meanings of r 23 , r 31 , r 10 , r 20 , and r 30 are 
similarly defined. Using the underlined r instead of pure r is simply for notational purpose; 
we will in a moment transform the coordinates into a more convenient form. Under the 
current notation Vq is given by 

Id d Id 2 

V-l = ^ (£o-f ) + ^4 (50) 
r o or o dr r 5 d#g 

and V 2 , V 2 , and V 3 are understood similarly. 



Performing the dimensionless unit transform same as made before in Sec. Ill] , the 
Hamiltonian is now 

, 1 1 1 1 1 1 \ 
+(— + — + • 51 

(V» iy» ty tf* V* V 

L\2 L.23 1-31 LlO L-20 L-30 

We introduce a new coordinate system 



ri 
r 2 



E 


+ Ei 


+ E 2 


+ E 3 


_ Eo 


+ Ei 


4 

-E 2 


-E 3 


_ Eo 


-Ei 


2 

+ E 2 


-Es 


_ E 


-Ei 


2 

-E 2 


+ E 3 






2 





rs = ^°~ M :^"^ 3 , (52) 



which are treated as vector transforms. 

These transforms result into the following relations, 



r 


- r x 


= E 2 + Es, 




r 


- r 2 


= E3 + E1, 




r 


- r 3 


= Ei + r 2 , 


(53) 



and 



r 3 -r 1= I3-T!. (54) 

With the transformed coordinate system in place, the Hamiltonian now takes a new 
form, 

h = -(ivj| + v? + vl + v3) 

/ 1 1 1 1 1 1 ^ 

+ — + — + ). (55) 

ri2 r 23 r 31 r 1+2 r 2+3 r 3+1 



where for example 



Id d Id 2 
V o = -^-( r o^-) + -^2, (56) 



and 



^12 = V r i + r i - 2r 1 r 2 cos6 l i2, (57) 



while ri + 2 is the short-hand for |ri + rgj or Jrf + r\ + 2rir2 cos^i2- 

ro now refers to the center-of-mass (cm.). After discarding the irrelevant first term on 
RHS of Eq. (|55|), we reduce the internal Hamiltonian in the cm. frame to 

/111 1 1 lx 

+ (— + — + ), (58) 

ri2 r 23 r 31 r 1+2 r 2+3 r 3+1 

which is effectively a three-body problem. 

B. Spin-singlet trial wavefunction of Is, 2p + , and 2p states 

It is easy to see that if the trial wavefunction for the three-body problem is antisym- 
metrized, the corresponding wavefunction for the our original four-body problem will be 
antisymmetric as well under the permutation of the original electron coordinates. Now we 
want to see how to make a three-body antisymmetrized wavefunctions. 
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The trial wavefunctions will be initially made only from the first four eigenstates for the 
purpose of simplicity 

[2 

| Is) = W— A exp(— Ar), 



|2p±) = V2^ A2rexp( ~Y )exp(± ^ ) - 



(59) 



which are copied from the previous section for quick reference here. As a notation conve- 
nience which will be extensively used later, we also define 



l 2 "> = V2^ A2rexp (-T ) ' (60) 

Note that (2p\2p) = 1 although \2p) is not necessarily orthogonal to other eigenstates. 

We now try to make an antisymmetrized wavefunction \1/ out of Eqs. \& should 

satisfy our main object, (^\H\ty) / {^\^) < 0. 

The nice symmetry property of the Hamiltonian helps us rule out \I> only made of s 
states. To see this, we need to understand that the kinetic terms in Hamiltonian are positive 
definite, or (^| - (Vf + V 2 + V"§) |^> /<^|^> > 0. Thus the only possibility left for a negative 

energy state is by a negative Coulomb term (\I>|(— + — + - — )|*)/(^|^). 

We only need to consider terms like — ^j)l^ r )/(^ r l^ r ) because the total Coulomb 

energy is just three times as much because of the symmetry of \& in r 1; r 2 , r 3 . Now we show 
and exploit a further symmetry hidden in r\ 2 and ri +2 . 

Note that 



T12 — yrf + r| — 2rir 2 cos i2 , (61) 



r i+2 = yr\ + r\ + 2r 1 r 2 cos#i 2 (62) 



'r\ + r\ — 2rir 2 cos(7r + 6> 12 ). (63) 
Both T\2 and ri +2 are periodic in Q\ with a period 27r. By the same argument we made 



in Sec. [TV] for two-body Coulomb integrals, we can integrate with respect to Q\ 2 and 9 2 , 
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both from to 2n. However r i+2 is identical to r\ 2 except for a phase shift of is in 6\ 2 . Since 
cos6>i2 has a period of 2i\ in 9 12 , (*|^|^)/(*|*) would be equal to (*|^|^)/(*|*) if 
\1/ does not have any angular dependency such as those given by pure s states. The total 
Coulomb energy would be and we have no way whatsoever to obtain a negative energy. 
We will give a more general theorem later and exploit similar properties to greatly simplify 
our calculation. 

So one of the states chosen in \I> has to be \2p + ) or \2p~). From our trial-and-error 
experience, we find that it is generally better to have two 2p states and one s state. 

As the first trial we want to test wavefunctions made out of |ls), \2p + ) and \2p~). 

By using variational approach, we also want to investigate the spin structure of such 
3E+1H complex. Since there is no spin-orbit coupling term in the Hamiltonian, we need 
only to specify the total spin of three electrons for any trial wavefunction. Basically there 
are two possible options for three electrons: either they are all lined up in the same direction 
or one of them is against the other two. The hole is also a fermion, whose spin tends to be 
opposite to at least one of the electrons. Thus the total spin of this complex is either or 
h, corresponding to spin singlet or triplet states. We would consider this explicitly in our 
calculation. 

First let us start with three orthogonal states \ls \), \2p + f) and \2p~~ j). The arrows 
designate spin states. The complex is now in a spin singlet state because the mixed spin 
wavefunction | |) = -^=| 1") + -^=| J,) should combine with that of the hole to give a singlet 
for one electron- hole pair. Writing | J) also means that we do not care about the exact spin 
state; it can be either up or down. 

The corresponding Slater determinant is 

|lsiti) |2p+ti) |2prli> 
\jl \ls 2 \ 2 ) \2p+h) |2p 2 -| 2 ) 
I lasts) \2pth) |2p 3 " Is) 
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\2pt h) |2p 3 Is) 



\2 P f ti) |2 P r ii) 
\2 P i h) 12P2 12) 



|l*s Is) 



|2 P + ti) |2pr ii> 

|2^ Ta) \2ps Is) 



= Ji{|2p 2 2p 3 ){exp[z(e 2 - 3 )]| T 2 | 3 ) - exp[^ 3 - 2 )]| T 3 | 2 )}|lsi |i) 

+ \2p 1 2p 2 ){exp[i{6 1 -0 2 )]\ T1I2) -exp[z(^ 2 -^)]| T 2 |i)}|ls 3 | 3 ) 

- |2 Pl 2p 3 ){exp[i(0i ~ 0*)]\ Tils) - exp[z(0 3 - 9 1 )}\ Ui)}\^2 I2)}, 

where for example |i indicates a mixed spin state for electron 1 while 2p2~ says that electron 
2 is in a 2p + state and so on so forth. Subscripts 1, 2, and 3 as usual are electron indices. 

There will be 6 terms in this antisymmetrized trial wavefunction, which would render 
follow-on manipulation difficult. Here we do a small trick by adding to the above the 
following Slater determinant, which is similar to Eq. fl6"4|) except for the switch between 2p + 
and 2p~ states. 

li> |2j?r to Mil) 

|ls 2 I2) |2p 2 T2) \2pt | 2 ) 
|ls 3 Is) |2ps h) \2pt Is) 

|2p 2 T2) |2p+ | 2 ) 
|2p 3 Ts) |2p+ | 3 ) 



(64) 



Ii> 



|2^r ti) |2rf ii> 

|2P2 h) \2pt I2) 



|l*s Is) 



\2Pi Ti) |2rf |i) 
\2P3 Ts) |2p+ | 3 ) 



|l*a I2) 



~{|2p 2 2p 3 ){exp[i(0 3 - 9 2 )]\ t 2 | 3 ) - exp[z(fl 2 - 3 )]| Tsl 2 )}|lsi |i) 
+ |2 Pl 2p 2 ){exp[z(0 2 -0i)]| Til 2 }-exp[i(0i-0 2 )]l T 2 li>}|ls 3 Is) 



- |2 Pl 2p 3 ){exp[i(0 3 ~ Tx| 3 ) - expfifa - 3 )]| Tali)}!^ h)}- 
The result is the following after normalization, 

tt = -^{12^2^133)008^21^3) 
+|2p 2 2p 3 ls 1 )cos0 23 |<S 2 \ 1 ) 

+ |2p32pil S 2)cOS^l|5 3 a li2 )}, 



(65) 



(66) 



22 



where the spin wavefunction |«S° 2 ,3) = "^(1 T1I2) _ I T2 -t 1 ) ) I $3) which is antisymmetric in 
the first two indices 1 and 2, and note that 

(«^12,3 1 ^23,1) = (^12,31^31,2) = (^23,ll^31,2) = ~2' 

is already normalized and thus (^\^) = 1. The kinetic energy of this trial wavefunction 
(^| - V 2 - V 2 - is easy to work out 

<*|-V 2 -V 2 -V 2 |*) = yA 2 . (68) 

Actually this can be figured out without explicit calculation. We have three particles in Is, 
2p + and 2p~ states whose kinetic energy are respectively A 2 , A 2 /9 and A 2 /9, thus the total 
kinetic energy should be the sum, which gives llA 2 /9. 

As for the Coulomb energy, we need only to calculate this term 

(*|- -|*>, (69) 

because 

(*|— " — 1*> = <*| — " — 1*> = (*|— " — 1*>, (70) 
ru r 1+2 r 23 r 2+3 r 31 r 3+1 

considering symmetry in electron permutation. 
We then write explicitly 

<*|— -— 1*> 
n 2 r 1+2 

= U(2 Pl 2p 2 ls 3 \ cos 2 9 12 (- — )|2 Pl 2p 2 ls 3 ) 

6 7-12 r 1+2 

+ (2p 2 2p 3 l Sl \ cos 2 23 (- —)\2p 2 2p 3 l Sl ) 

ri2 r 1+2 

+ (2p 3 2 Pl ls 2 \ cos 2 3 i(- — )|2p 3 2 Pl l S2 ) 

ri2 r 1+2 

- (2p 2 2p 3 lsi\ cos6»i 2 cos6» 23 ( )|2pi2p 2 ls 3 ) 

n 2 r 1+2 

- (2p 3 2pils 2 | cos 6»i2 cos6» 31 ( )|2pi2p 2 ls 3 ) 

rn r 1+2 

- (2p 3 2 Pl ls 2 \cos9 23 cos6 31 (- —)\2p 3 2p 2 l Sl )}. (71) 

n 2 r 1+2 
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For obvious identities from symmetry 



(2p 2 2p 3 l Sl \cos 2 6 23 {- l —)\2p 2 2p^ Sl ) 

ri 2 r 1+2 

= (2p 3 2 Pl ls 2 \cos 2 e 31 (- — )|2p 3 2 Pl Ls 2 >, 

ri 2 r 1+2 

(2p 2 2p 3 lsi| cos 9 12 cos 9 23 ( )\2p 1 2p 2 ls 3 ) 

= (2p 3 2p 1 ls 2 \ cos 6>i2 cos 6>3i ( )|2pi2p 2 ls 3 ), 

ri 2 r 1+2 

Eq. (|7l"D is simplified to 

— 1*> 

ri2 r 1+2 

= -{(2 Pl 2p 2 ls 3 \cos 2 9 12 (- —)\2 Pl 2p 2 ls 3 ) 

t> r i2 r 1+2 

+ 2(2p 2 2p 3 lsi| cos 2 9 23 {- — )|2p 2 2p 3 lsi> 

ri2 ri +2 

- 2(2j9 2 2p 3 lsi| cos6» 12 cos6'2 3 ( )|2pi2p 2 ls 3 ) 

^12 r 1+2 

- (2p 3 2pils 2 | cos 023 cos 31 (- — )|2p 2 2p 3 lsi>}. (72) 

ri2 t*i +2 

The first term on RHS of Eq. (|72|) will be zero. To understand this, we prove the 
following theorem. 

Theorem: Suppose that /i(Vi, r 2 , r 3 ) and fz(ri, r 2 , r 3 ) are functions only of r 1; r 2 , and r 3 
(z.e., do not depend on 9i, 9 2 , and 6 3 ). Also assume that (/iJ^I/2) exists where fi and f 2 
are short-hands for fi{ri,r 2 ,r 3 ) and f2{fi,r 2 ,r 3 ). Then we have 

(Al COS" U (- — )l/2> 

^12 r 1+2 



0, if n is even; 

(/ 1 |2cos^ 12 (i)|/2),ifnis odd. 



(73) 



Proof. 



(/l|00S n M— -— )l/2> 
^12 ^1+2 



/"OO /"CO /•go /-27T /-27T /-27T j j 

/ ndn / r 2 dr 2 r 3 dr 3 dOj d6 2 ^3/1/2 cos" M ) 

JO JO JO JO JO JO 7*12 ?"l+2 

/ ri drx r 2 dr 2 r 3 dr 3 d6 3 d9 2 d9 12 f 1 f 2 cos n B l2 { ) (74) 

jo jo Jo Jo Jo Jo r\2 ri +2 
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where in the second step we have changed variables from 6\ to 0i 2 for the same reason that 
we discussed before in Sec. fVj. We consider first the integration with respect to 0i 2 in Eq. 
([T^), or the following expression 

2tT 11 

de 12 cos n e 12 ( ) 

ri2 ri +2 

2tt 1 1 

de 12 cos n e 12 ( -, , )■ (75) 



ijr± + r 2 - 2r 1 r 2 cos9i 2 yj r\ + r\ - 1r x r 2 cos(7r + i2 , 

The last term on RHS of Eq. ( fT5|) is 

^0i 2 cos n 0i 2 ; 1 (76) 

rf + r 2 - 2r\r 2 cos(7r + i2 ) 



Changing variable from 12 to an arbitrary x = 7r + 12 in this integral, we have 

2tt 1 
(2012 COS™ 012- 







r\ + r\ - 2r x r 2 cos(7r + 0i 2 ) 



3tt 1 

dx cos™ fx — it) 



r{ + r 2 — 2r\r 2 cos x 

' W dx(-l) n cos n x - / 1 (77) 

r i + r i — 2rir 2 cos x 



Because the integrand on RHS of Eq. ( 77 ) is a function of x with a period 2ir, integrating 
from to 7r is equal to integrating from 27T to 37T, or 

r37r 1 



riir 

/ dx{-l) n cos n x 

J2-K 



rf + r 2 — 2rir2 cos a: 



dx{-l) n cos" x 1 (78) 

*' w 7*1 + r| — 2rir 2 cos x 



or alternatively 



3tt 1 

dx(-l) n cos n x 



r i + r l — 2rir 2 cosx 

2?r 1 



- / * dz(-l) n cos n x (79) 
J0 Jrf + r 2 — 2rir 2 cos a; 



After Eq. is substituted in Eq. ([77]) 
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2tt I 

d9 12 cos n 9 12 - 



r l + r 2 — 2r l r 2 C0S(7T + 9 12 ) 



2tt 1 

dx(-l) n cos n x- 



ii t l r 2 _|_ r | _ 2r\T2 cos X 

2tt I 



- d9 12 (-l) n cos n 9 12 (80) 

7,1 ' ' t\ — (— 9™2 — lT\T 2 COS #12 



where in the second step variable x is relabelled as #12. 



Inserting Eq. ( p0|) into (|75|) , we have 



27T 11 

d012 COS" M 



ri 2 ri +2 

2tt 1 (-1)" 

d# 12 cos n # 12 [ 



" yr 2 + rf — 2r x r 2 cos 6^2 yr 2 + rf — 2rxr 2 cos#i 2 

0, if n is even; 

/o 2 "^ 12 ^^,ifnis odd. 



(81) 



Then we substitute Eq. (^l|) back to Eq. (|74|) to give the theorem. This concludes our 



proof. 

The second term on RHS of Eq. (|72D or 2(2p 2 2p 3 lsi| cos^ 3 (^ - ^)\2p 2 2p^lsi) has 
the factor 

cos 2 #23 = [cos 9 2 cos 6*3 + sin 9 2 sin 6*3] 2 
= cos 2 6*2 cos 2 6*3 + sin 2 9 2 sin 2 # 3 

+2 cos 6*2 cos #3 sin 9 2 sin # 3 , (82) 

where last term on the right gives after integration with respect to #3, while the first two 
terms, combined with ^ — — - are after integration with respect to #12 from our theorem, 
because there is no odd power of cos #12 in the integrand. 

The third term on RHS of Eq. (|72|) is after integration with respect to #3 and we have 
only the last term left. 

Because 

COS #23 COS #31 
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[cos 6*2 cos 83 + sin # 2 sin 63] [cos 6*3 cos 6\ + sin 63 sin #1] 
cos #2 cos 9\ cos 2 #3 + sin 82 sin #1 sin 2 O3 

+ cos 6*3 cos 6*i sin 6* 2 sin # 3 + cos 62 cos #3 sin 63 sin #1 , (83) 



which after integration with respect to 83 gives 



7r(cos 62 cos 9 1 + sin 62 sin 8\) 

= 7TCOS#i 2 . (84) 



Thus 



■ (2p 3 2pi ls 2 1 cos #23 cos 3 i (— — ) 1 2p 3 2p 2 lsi ) 

ri2 n+2 

: _( 2pi l S2 |£^|2p 2 l Sl ) (85) 

fl2 



which is worked out in the same fashion as integral (|30[) in section |IV| to give 



(2p 1 lg 2 | CQS ^ 12 |2p 2 lgi) = -0.0172573A. (86) 
f 12 



The total Coulomb energy is thus — 0.0345146A and we have 



E 



11A 2 



- 0.0345146A 



9 

= —(A - 0.0141196) 2 - 0.00019936, (87) 

from which the upper bound for ground-state energy is therefore —0.00019936. 

This confirms our conjecture that 3E+1H indeed forms a quasi-bound state. Here quasi- 
bound merely means that the 4 particles are relatively localized or near one to another and 
the total energy of this system is lower than when they are infinitely apart. Note that the 
4 particle as a whole could still move freely with an effective mass of 4m. 

We proceed on to try other wavefunctions and want to improve our result. 
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C. Spin-triplet trial wavefunction of Is, 2p + , and 2p states 



The next logic choice is still Is, 2p + , and 2p states but with a spin-triplet configuration, 

vi 1 



2 

# = - 7 ={\2p 1 2p 2 ls 3 )sm9 12 \S s 123 ) 



+ \2p 2 2p 3 ls 1 ) sin 6>2 3 |«S| 3! i) 

+ |2p 3 2pils 2 ) sin#3i|S 3 s li2 >}, (88) 

where the normalized spin wavefunction |iSf 23 ) is symmetrized in indices 1 and 2, and 
describes a spin triplet state for electrons 1 and 2. |tSf 23 ) may be written | ti?2 J3) — 2 I T1T2 
)(l T3) + I I3)) but can also take other forms, for example | I1I2I3) or |^(Til2 + T2I1) U)- 
Three forms of \Sf 2 3 ) all satisfy (£12,3 1^12,3) = 1 an d (^12,31^23,1) = 1/2 and so on. 

This trial wavefunction gives exactly the same result as the spin singlet trial wavefunction, 



E 



11A 2 
9 



0.0345146A 



= y (A - 0.0141196) 2 - 0.00019936, (89) 

which gives E min = —0.00019936. We are not yet able to tell whether the spin singlet or 
the triplet state is more stable. 



D. Spin-singlet trial wavefunction of modified Is, 2p + , and 2p states 

We know in such a 3E+1H complex, the strong repulsive Coulomb force is a huge penalty 
we have pay attention to, thus it would be nice if our trial wavefunction is small when two 
electron coordinates are close. Take the spin-singlet wavefunction for example, 

# = -^{|2 Pl 2p 2 ls3> cos 9 12 \S? 2>3 ) 

COS #23 1 ^23 1 ) 

+ \2p 3 2 Pl ls 2 )cos9 31 \S«)}, (90) 
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if we can modify \I/ so that it is small when two electrons are close, the total energy then 
might be lower. So we choose the following, 

= -^={|2 Pl 2p 2 l S 3>(l - cos £ 12 ) 1^2,3) 
+]2p 2 2p 3 ls 1 )(l-cos^3)|<S 2 a 3 jl ) 

+ |2p 3 2p 1 l S2 )(l-cos^i)|5 3 a li2 )}, (91) 

where the first term on RHS will be zero when i2 = 0, or 6\ = 2 , or in other words, when 
electrons 1 and 2 are near. 

Note that after this modification, the three particles are no longer in pure Is, 2p + , and 
2p" states, but we may still use these words for notational purpose, and call them modified 
Is, 2p + , and 2p~ states. 

Then 

(*|*) =2(2 Pl 2p 2 ls 3 |l -2cos# 12 + cos 2 9 12 \2 Pl 2p 2 ls 3 ) 

-2(2pi2p 2 ls 3 |l - cos 0i2 - cos0 23 + cos6>i 2 cos0 23 |2p 2 2p 3 lsi) 
= 3-2{(2p 1 ls 3 |2p 3 lsi) + (2pi2p 2 ls 3 |cos0i 2 cos0 2 3|2p 2 2p 3 lsi)}, (92) 

where we have used identities 

(2p 1 2p 2 ls 3 \l -2 cos 0i2 + cos 2 0i2|2pi2p 2 ls 3 ) 
= (2p 2 2p 3 ls 1 \l - 2 cos 23 + cos 2 02 3 |2p 2 2p 3 lsi) 

= (2p 3 2pils 2 |l -2cos0 3 i + cos 2 03i|2p 3 2pils 2 ), (93) 

and 

(2pi2p 2 ls 3 |l - cos 0i2 - cos 23 + cos 0i2 cos 23 1 2p 2 2p 3 1 Si) 

= (2p 2 2p 3 lsi\l - cos 23 - COS 3 i + cos0 23 cos03i|2p 3 2pils 2 ), (94) 

and the like. 

The exact value of (^|^) is found to be 90/32. 
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Next we want to find out the kinetic energy 

(*| - V? - v 2 - v 2 |*} 



where 



Because 



we have 



3(*|- V?|*) 

3(Vi*|Vi*), (95) 



d 1 d 

Vi = ri — + 0^ — . (96) 



Vi|*> = -={(- - -)|2 Pl 2p 2 ls 3 >(l - cosM^, 



-A|2p 2 2 P 3l Sl )(l-cos^3)5 2 a 31 



+ - ^)|2p 3 2j?ils a >(l -cos0 3 i)5?i, 2 }fi 

+ -4 — {\2 Pl 2p 2 ls 3 ) sm9 12 S^ 3 
V 6 n 

- \2p 3 2 Pl ls 2 ) sin^i<S 3 a li2 }, (97) 



(Vi*|Vi*> 

4 1 2A A 2 

= {2(2 Pl 2p 2 |(- - — + — )(1 - 2cos# 12 + cos 2 ^ 12 )|2 Pl 2p 2 ) 
6 rf 3ri 9 

+ (2p 2 2p 3 |A 2 (l + cos 2 23 - 2cosfl 23 )|2p 2 2p 3 > 

-2(2 J ? 1 ls 3 |(^--)|ls 1 2p3> 
3 rx 



1 2A A 2 

(2 Pl 2p 2 l S3 |(- - — + — )|2 Pl 2p 3 ls2> 
rf 3ri 9 



■ sin 2 0i2 1 , , , 

+ 2 (2 Pl 2p 2 1 ^ 1 2 Pl 2p 2 > }, (98) 

where we have used symmetry identities similar to Eqs. ( |9"3"| ) and (|9~3|)to merge some terms. 
Those terms which go to after integration of one variable, say 6 2 , are dropped as well. 
After somewhat tedious but straightforward calculation, the total kinetic energy is ^r^A 2 
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The Coulomb energy of the trial wavefunction is calculated as usual though with added 
complexity. We just give the result here, — 0.566523A. 
Thus 

(V\H\W) 



E 



1.15638A 2 - 0.201430A 



= 1.15638(A - 0.0870951) 2 - 0.00877178, (99) 

which gives E min = —0.00877178, a much improved result over previous cases. 

Also recall that we have spin-singlet state in this case. In contrast, the spin-triplet trial 
wavefunction of Eq. (|90|) cannot be so modified as Eq. (f)TD because the first term on RHS 
of Eq. ( |9~0"D is already small when electrons 1 and 2 are close. This suggests that antiparallel 
spins be energetically more favorable than parallel ones. We will from now on focus our 
attention to spin-singlet state only. 

E. Spin-singlet trial wavefunction of modified 2s, 2p + , and 2p states 

The next trial wavefunction and the final one in this paper is made of modified 2s, 2p + , 
and 2p~ states, 

* = ^={\2 Pl 2p 2 2s 3 )(l - cos £ 12 ) 1^,3) 
+ |2p 2 2p32si)(l-cos^3)|5 2 a 3il ) 

+ \2p 3 2 Pl 2s 2 )(l -cos6 31 )\S^ 2 )}. (100) 

The treating of such \1/ is similar to the previous case, and we just list the final results 
as follows, 

(mm 



E 



0.632997A 2 - 0.168735A 



= 0.632997(A - 0.133283) 2 - 0.0112447. (101) 
Then the upper bound for ground-state energy is now improved to —0.0112447. 
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VI. STABILITY AGAINST ELECTRON-HOLE RECOMBINATION 



The stability of a trion or a tetron against electron-hole recombination or annihilation ac- 
tually is not a concern if the electrons and holes are from semi-metallic band-edge. Whether 
electrons (e) and holes (h) should recombine or not depends on the band structure. Here we 
just give a very simple example: semimetals. Electrons and holes coexist peacefully without 
annihilating one another, in semimetalsEl such as Bi, As, Sb. The reason is as follows and 
shown in Fig. 4. 




FIG. 4. Mechanism for the stability of electron-hole pairs, (a) Typical band structure for semi- 
conductors. Electron states have higher energy than hole states. As a result, electrons recombine 
with holes in a short time. This reasoning also applies to metal provided the relevant conduction 
band is partly filled, (b) and (c) Band structures of semimetals. Electron states have lower energy 
than hole states and the recombination of electrons and holes is unfavored in energy. In all figures, 
Ep is the Fermi energy. 

In typical semiconductors, the conduction and valence bands are separated by a positive 
energy gap as shown in Fig. 4(a). By jumping from the conduction to the valence band which 
is at a lower energy (i.e., recombining with a hole), the thermally excited electron lowers 
its energy. This process thus occurs spontaneously to the e-h pair and makes it unstable. 
(The same process can be interpreted as a hole jumping up and note that in contrast to 
electron, the higher is the hole in the band diagram, the lower is its energy.) For semimetals 
with different band configurations, as typified in Fig. 4(b) and (c), the situation is however 
totally different. If the electron were to jump to the hole state and then recombine, as 
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indicated by the arrows, it would have to take extra energy to reach there. This process 
is clearly unfavorable in energy and such a mechanism secures the stability of electron-hole 
liquid against recombination in semimetals. Thus we expect the proposed tetrons to find 
their best playing ground in solids with semi-metallic band structures. 

VII. STABILITY AGAINST DISSOCIATION 

As mentioned before a tetron is not stable against the dissociation of one of its electrons 
into infinity in vacuum. To understand this, let us perform a thought experiment. The 
theorem by Lieb states that the ground-state energy of a 3E+1H complex must be higher 
than that of a trion plus an electron which is assumed to be infinitely far away, having 
arbitrarily small kinetic energy. Thus after we place a quasi-bound state of 3E+1H in 
vacuum, the complex tends to lower its energy by moving one of the electrons away to 
minimize the energy. This corresponds to the dissociation (detachment) of one electron as 
shown in Fig. 5. 

° 3 electrons + 1 hole complex 
(tetron) in vacuum 

Dissociation into a trion and an electron 

This electron 
moves to infinity 
o o 

FIG. 5. The reason why tetron (3E+1H) is not stable in vacuum. 

However if the tetron is not in vacuum, for example, if many other electrons and holes 
are present nearby, then there is hardly any gain in energy for the whole system even if an 
electron dissociates from one tetron; this electron would simply increase Coulomb energy 
with other electrons anyway even if it decreases the energy of the original tetron as shown in 
Fig. 6. Such cases may happen to condensed matter systems where charge carrier densities 



o 
o 

V 
o 

Trion o 
o 
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are governed by electronic structure, or some types of external controls such as dopings. The 
electrons cannot leave the holes for infinity because statistically there must be some electrons 
close to the holes. When the densities of electrons and holes are low, the situation is more like 
a vacuum and therefore tetrons are not stable. With an increasing charge density, however, 
the story could be increasingly different. Then if the densities exceed a certain threshold, 
the whole system would lower its energy by taking advantage of the favorable geometric 
configuration offered by tetrons, provided the electrons outnumber holes by a ratio near 
or higher than 3. Then the tetron would be more stable in this high density electron-hole 
liquid. Actually Nature may have already suggested a similar mechanism for us: doubly 
charged negative ions until recently have not been seen in gas phase, or in isolation. In 
contrast, they have long been known to exist in condensed matter systems such as liquids 
and solutions^!. 

If the densities are to further increase by too much, the energy of the tetron state will 
no longer be negative after a second threshold. The reason is simple, the gain in Coulomb 
energy increases as a -1 while the penalty paid for the kinetic energy scales as a~ 2 , where 
a is the average distance between particles. Then we would see a system close to ordinary 
metal because the quasi-bound state is no longer possible. 

Here we just give the above qualitative argument to show how tetron states might be 
materialized. Further studies are needed in a more quantitative and rigorous way. 
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Where to dissociate? The electron cannot move to infinity. 



FIG. 6. The reason why tetron (3E+1H) could be stable in the electron-hole liquid system. 



VIII. POSSIBLE CONNECTION TO HIGH- Tc PAIRING MECHANISM? 

The tetron state acts in a whole as a boson because it starts with 4 fermions. This 
is significant because the net charge is — 2e (2e if we start with 3h+le which is a more 
realistic scenario for most copper-oxide superconductors), and we are now left with a Cooper 
pair in effect. Moreover our current variational result implies a singlet state for the spin 
part of wavefunction, which is exactly what has been found in the Cooper pairs of high-T c 
superconductors!!. This entices us into the possible connection between the tetron state 
and high-T c superconductivity, although this model itself could be oversimplified without 
considering the crystal structure of the host material. 

Charged fermion pairing models based on exchange of excitons have been proposed before 
by some notable theoristsSill, often in a paradigm of perturbative approaches. What we 
have developed here could be a non-perturbative way of treating similar problems. 

High-T c superconductivity found in copper oxides or cuprates are likely from the Bose- 
Einstein condensation of pre-formed pairsS In this regard, the tetron complex reasonably 
fits into such a picture, with its small size comparable to the coherence length of cuprate 
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Cooper pairs. Also high-T c is intimately related to low dimensionality. Recall that the 
ground-state energy of a 2D exciton is four times as low as that in 3D, while T c decreases in 
cuprates when the superconductors are increasingly doped from a 2D conductor to 3D one. 

There are mounting evidences showing the coexistence of electrons and holes in cuprate 
superconductors from key experiments such as angle-resolved photoemission spectroscopy^ 
and transport properties@S. Finally the qualitative discussion given in Sec. [VII] actually 
quite fits to the doping- dependent phase diagram of cuprate superconductors from under- 
doping to overdoping. Although the case at this moment is far from certain, the author 
sincerely hopes that this paper may generate some interests along this line of reasoning. 



IX. CONCLUSION 

In conclusion, various trial functions for 3E+1H complex (tetron) were worked out by 
rigorous variational method to give an upper limit on the ground-state energy. The calcula- 
tion suggests quasi-bound states for such a complex. The possible connection of such states 
to a pairing mechanism of high-T c is discussed. 
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